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Abstract 

The influence of the HBT effect on the multiplicity distribution 
and charge ratios of independently produced pions is studied. It is 
shown that, for a wide class of models, there is a critical point, where 
the average number of pions becomes very large and the multiplicity 
distibution becomes very broad. In this regime unusual charge ratios 
(" centauros" , " anticentauros" ) are strongly enhanced. The prospects 
for reaching this regime are discussed. 

1 Introduction 

It is now well established that the HBT correlations |l[ influence significantly 
the momentum distribution of particles created in high-energy collisions. The 
effect on multiplicity distributions and on particle ratios, however, although 
predicted theoretically 0- ||, has not yet been found in accelerator experi- 
ments ||. In view of the increasing interest in measurements of multiplicity 
distributions in future collider experiments, we found it worthwhile to reex- 
amine this problem once more. 

*Also at the Nuclear Physics Institute, Krakow, Poland 

t Address: Reymonta 4, 30-059 Krakow, Poland; e-mail: bialas@thp4.if.uj.edu.pl, 
zalewski@chall.ifj.edu.pl 



1 



We begin by reviewing the basic ideas of the HBT effect, as applied to 
processes of particle production. This will permit to explain our assumptions 
and to introduce the notation[|. 

Let i/; (qi,q2, ■■■q n ,a) = ipo(q, a) be the probability amplitude for the 
production of n particles with momenta [qi, q 2 , q n ] = q calculated ignoring 
the identity of particles. Here a denotes a collection of all other quantum 
numbers which may be relevant to the process in question ( e.g., the momenta 
of other particles which we do not wish to consider explicitly in a "semi- 
inclusive" measurement). The density matrix 

Po(Q, Q') = J daip (q, a)^l{q' ', a) (1) 

gives all the available information about the system in question. In particular, 
the momentum spectrum of particles is 

^o(q) = J da | ip (q,a) \ 2 = p (q,q). (2) 

In the following we shall assume that flo(q) is normalized to 1: 

J dqn (q) = 1. (3) 

Suppose now that the particles are identical. In this case the states with 
different permutation of particle momenta are non-distinguishable and the 
wave function is a sum over all permutations 

= J2^o(qp,a), (4) 
p 

where qp is the set of momenta [qi,q<2, ■■■,q n ) ordered according to the per- 
mutation P of [1, 2, n]. Using (|l])-(§) we obtain for the distribution of 
momenta of identical particles 

If 1 

fi (<?) = — / da | ip(q, a) \ 2 = — ^ p (q P , q P ,). (5) 

Tl. J Tl. p pi 

The factor — x takes care of the fact that the phase-space for n identical 
particles is n! times smaller than the phase-space of the non-identical ones. 

1 We will follow the approach of Q. 



2 



Eq.([5]) summarizes the effect of the identity of particles on the observed 
spectra. It is seen that to evaluate this effect it is not enough to know the 
spectrum fto(q). The full density matrix po(q,q') is necessary. Inverting 
this statement we observe that the measurements of the spectra of identi- 
cal particles provide information on the density matrix po(q, q'), inaccessible 
otherwise. This is, in fact, the reason why these measurements are are so 
attractive. 

There are three points worth observing in 

(i) When momenta of all particles are equal, we obtain the simple result 

&{qi = <?2 = ••• = q n ) = n!O (<?i = q 2 = ... = q n )- (6) 

(ii) The normalization of fl(q) is different from that of £lo(q). Generally 
we have 

/ dqn(q) > J dqn (q) (7) 

Since the relation between fl(q) and fi (g) depends on n, the HBT effect 
affects not only the shape of momentum spectrum but also the multiplic- 
ity distribution. This aspect of the problem is the subject of the present 
investigation. 

(hi) If particles are emitted in a pure state, i.e., if Po(q,q') = V'cKQO'^oty) 
and if ipo(q) is symmetric with respect to interchange of any pair of momenta, 
we have 

Q(q) = n\tt (q) (8) 

a really dramatic result .0 

2 HBT phenomenon in uncorrelated emission 

Consider now a system of n particles emitted independently. If we ignore the 
identity of particles, independent emission implies that the density matrix 
factorizes 

n 

po{q,q') = Y[po{qi,q'i)- (9) 

i=i 

2 We emphasize again that flo(q) is the spectrum calculated with identity of particles 
being ignored. Some authors ||, while discussing the HBT effect for pure states, include 
the identity of particles already in calculation of Qq(q). In this case one obtains of course 
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Introducing this into ([5]) we have 

1 n 

"(?) = rrEIlA)((?p)i ) (?p')0- (io) 

Here we are interested in the integral 

W n = [ dqn(q) = E / ft {d\ lPo (q u (g P )0) . (11) 
J p J i=i 

To calculate W n we observe that, for each permutation P the integral on 
the right hand side of ([□]) factorizes into a product of contributions from all 
the cycles of P (as is well-known, each permutation can be decomposed into 
cycles). Let us denote the contribution from a cycle of length k by We 
have 

C k = J d 3 q 1 ...d 3 q k p (q 1 , g 2 )Po(<?2, qa)-Po(Qk-i, Qk)Po(<lk, ?i) (12) 

It follows from ([3]) that C\ — 1 . For k > 1, C k depends on the form of 
Po(q,q') and cannot be calculated without further assumptions. One can 
prove, however, that all Ck are positive. Indeed, one sees from ( Jl2| ) that 

C k = Tr[ Po } k . (13) 

Since po(q, q'), being a density matrix, has only non- negative eigenvalues and 
trace one, C k > 0. 

The rest of the calculation is just combinatorics. 

We observe first that any two permutations which have identical parti- 
tions into cycles give equal contributions. Denoting by n k the number of 
occurences of a cycle of lenght k in the set of permutations considered, the 
contribution from all of them can be written as 

( c k \ rik 

In the first equality the first factor is the integral, the second is the number 
of partitions of the n particles among the cycles, the third is the number of 
ways a cycle can be constructed from k particles and the last one corrects 
for the permutations of whole cycles. 
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W n is obtained by summing W' n only over permutations which have par- 
titions into cycles different from each other. This still leaves a large number 
of terms but -for large n- this number is much smaller than the original n\. 

Until now we have considered a fixed multiplicity. If the multiplicity 
distribution calculated with identity of particles ignored is given by Po(n), 
the correct multiplicity distribution of identical particles is 

= p o{ „)w a (15) 

For independent emission Po(n) is given by the Poisson distribution 

P Q (n) = e~ v — y (16) 
n! 

and we obtain an elegant formula for the generating function of the multi- 
plicity distribution: 

*(z) = £ P(n)z n = exp (± " V : lK7fc ) • (17) 

n \k=X K ) 

3 General discussion of multiplicity distribu- 
tions obtained from independent emission 

In this section we discuss the general properties of the multiplicity distribu- 
tions obtained from the Eq. (|TTD . 

First, we observe that using well-known properties of the generating func- 
tion we obtain from (|TT|) for the average multiplicity 

oo 

<n>=Y,v k C k , (18) 

k=l 

and for the correlation coefficients (cumulants) 

^It^' (19) 

All the cumulants are positive, because all the Ck are. This means in partic- 
ular that the distribution is always broader than the Poisson one. 
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Specific properties of the distributions defined by ( |TTD depend, of course, 
on the value of v and of the cycle integrals Ck- 

The first important example we would like to consider is when particles 
are emitted in a pure state, i.e. 

Po(q,q')=^(qW(q')- (20) 
It follows from ( |T2"D that C& = 1 for all k and the generating function becomes 

corresponding to the geometric distribution 

P{n) = (1 - v)v n , < n >= j^— (22) 

which, at the critical point v — > 1, exhibits the phenomenon of Einstein 
condensation. One sees from this example that the resulting multiplicity 
distribution has little to do with the original Poisson one, and the observed 
average multiplicity may dramatically differ from the initial v. As we shall 
see, this is a general phenomenon. 

Further discussion depends on the assumed shape of the single-particle 
density matrix. 

The evaluation of Ck is greatly simplified if one works in the basis where 
the density matrix po(q,q') is diagonal. 

Let us first discuss the case of a discrete eigenvalue spectrum. We have 

C* = XX Ci = £A m = l. (23) 

m m 

and thus the generating function of the multiplicity distribution can be rep- 
resented as a product 

Hz) = H* m {z) (24) 

m 

where 

1 - z\ m v 

are the generating functions of the geometric distribution (if the eigenvalue 
X m is degenerate, the corresponding factor in (j2iD appears g m times, where 
g m is the degeneration factor). 
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Thus we obtain the average multiplicity 



< n >= > > — (26) 

^ 1 - \ m u ~ 1 - X v 

and the cumulants 

X m v Y , , ... / An^ 



fr-U'E ifxr ^- 1 Hi^ (27) 



where Ao is the largest eigenvalue of po(q,q'). One sees that z/Ao = 1 is the 
critical point of the multiplicity distribution (c.f. @). 

Let us now discuss the case when the eigenvalue spectrum is continuous. 
Denoting by A the eigenvalues of po(q,q') and by cr(A) the spectral function 
we have 

C k = f a(X)X k dX. (28) 



The normalization condition (0) implies that 

d = [ a(X)XdX = 1. (29) 



One sees that the problem reduces to a discussion of a single non-negative 
function er(A) defined in the interval [0, 1] and satisfying the condition (]29]) . 
We also note that the eigenvalues equal to zero do not contribute to One 
can thus always add to cr(A) a term of the form a5(A) (with an arbitrary 
positive constant a) without changing the results. 

To be more specific, we consider the generic spectral function in the form 

for < A < Ao and zero otherwise. 
We obtain 

r _ , k ^ T(a + b + 3)T(a + k + l) 

Uk ~ A ° T(a + 2)T(a + b + k + 2) m 



and the cumulants 

^- r( p ( ::;^ ( ;^ ( :;y ^v fan+y+1;a+t+P+2; v).(32) 
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The hypergeometric function F becomes singular at uXq = 1 when p exceeds 
6+1. 

Three special cases are of interest: 

(i) b — > —1. In this case cr(A) — > — A ) and we have Ck = Aq 1 . 
Thus 

and we recognize the negative binomial distribution with the average < n >= 
1 _ u Xq1 , - For Ao — > we recover the Poisson distribution flTEp. As discussed in 
the previous section, Ao — > 1 corresponds to the pure state Cj. = 1. 

(ii) a = b = 0. Now we have = ^-Ao -1 an d one finds the following 
formula for the cumulants of the distribution 

K p = 2^-^^(-\-Y F(l,2;p + 2;A ^. (34) 
p+1 \l-\ Q vJ 

From the know properties of the hypergeometric function |J we deduce that 
at the critical point Ao^ — ► 1, < n > diverges logarithmically, whereas all 
other cumulants diverge as negative powers of (1 — Ao^). 

(hi) b > 0. In this case the average multiplicity approaches at the 
critical point Ao^ — > 1. The divergences appear only in cumulants of order 
p > b + 1. Such a distribution is characterized by a large probability of 
no particles at all, accompanied by a small probablity of production of very 
many particles. Schematically 

^(O)- 1 "^ P{n) = j^Pin) n>0. (35) 

Here N — > oo is the average of the distribution p(n). 

To summarize, the distribution becomes always singular when u\o (i.e. 
the product of the initial average multiplicity and the maximal eigenvalue of 
the density matrix) approaches 1. The character of the singularity at this 
critical point, however, depends crucially on the behaviour of the spectral 
function in the vicinity of Aq. 
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4 Gaussian density matrix 

In this section we consider the single particle density matrix of the gaussian 
form, discussed already in several papers by Pratt |@, [| f§ (see also ||) 



Po(?, q') = Px(q x , q' x )p y (q y , q' y )pz(qz, q' z ) (36) 



with 



.ii3Li__I p2/„-N2 



where 



P^^=[^ x ) e-^r-> RM , (37) 

q + = \{q + q')\ q' = q-q'- (38) 

Analogous formulae define p y and p z . As easily seen, is the average value 
of the square of the x-component of the particle momentum, and R%, is the 
average value of the square of the x-coordinate of the particle emission point. 
The uncertainty principle implies that for i — x,y, z, 

R*\ > \ (39) 

In order to determine the multiplicity distribution, we first find the eigen- 
values of the density matrix (|37|) . To this end we observe that the eigenfunc- 
tions of ( |3TD are of the formQ 

f m {q)=e^ 2 H m {^q) 1 (40) 

where H m (q) is the Hermite polynomial of order m. Using Q4"0| ) it is not 
difficult to find the eigenvalues: 

A m = A rst = A (l - \ x ) r (l - \ y ) s (l - X z )\ r,s,t = 0, 1, (41) 

where A = A z A y A 2 and for i = x,y,z 

x - = (TT-bo (42) 



3 This was pointed out to us by A.Staruszkiewicz. 
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Note that (^) implies that Ao < 1, as necessary. 

The generating function is given by ([24]), (pf) and the cumulants by (p7|) 
with A m given by (gTJ). 

Using (f4~l|) and following the arguments of the previous section we also 
obtain the elegant formula 

A^ 

° k = ¥ A " = [l-(l-A^][(l-l-A^][l-(l-A^]- (43) 
Together with (0), this gives for the cumulants 

which diverges, as expected, at vXq — > 1. 



5 Charge ratios 

It was first pointed out by Pratt || || [| that studies of the charge ratios may 
be an effective way to uncover the effects of HBT correlations in multiparticle 
systems. This issue can be readily treated using the methods developped in 
Sectons 3 and 4. 

We would like to discuss independent production of positive, negative 
and neutral pions. The main difficulty in formulating the problem is how to 
implement the constraint of charge conservation, as this clearly depends on 
the dynamics of the production process (a thorough discussion can be found 
in Here we restrict ourselves to two cases which, we believe, illustrate 
well the main point: the charge ratios obtained may drastically differ from 
those expected from the "uncorrected" distributions. 

In the first case the constraint of charge conservation is ignored altogether. 
This may be justified if the system of particles we consider is a small part of 
a very large system. For the generating function of multiplicity distribution 
we thus obtain simply 

z_, zq) = $(z + )$(z„Mz ), (45) 

where §{z) is is the generating function of the multiplicity distribution of 
one of the species. Introducing n c = n + + n_ we have 

$(z c ,z ) = $ 2 (z c )$(z ), (46) 
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From this equation one can obtain the full joint distribution of charged and 
neutral pions by the usual methods. For illustration we just quote the results 
for two extreme cases n = ("centauros") and n c = ("anticentauros"): 

P(n = 0) = $(0), P(n c = 0) = [$(0)f. (47) 

Using now the results of the previous section we have 

P(n = 0) = " Knv), P(n c = 0) = [P(n = 0)] 2 . (48) 

m 

One immediate consequence is that the production of "centauros" must be 
larger than that of "anticentauros". For the Gaussian density matrix we 
obtain, according to ([41|), 

P(n = 0) = J] [l - AoKl - A.) r (l - A,r(l - A,)*] . (49) 

rst 

The second case we consider is when the charged particles are produced 
in pairs, i.e., the distribution is of the form 

. \P(n c /2)] 2 P(n ) , . 

Thus for production of "centauros" we obtain the same formula as before. 
The generating function reads 

$(z C) z ) = *(z c )$(zo) (51) 

where 

T/ s En [P(n c /2)] 2 z^ 

V(z c ) = n ^ p 7 2 c ■ 52 
The probabilty of creating an "anticentuaro" is 

* (4=0) ™ >|4(1)f (53) 

and thus the production of " anticentauros" is enhanced as compared to the 
previous case. 

These results are illustrated in Fig.l, where the probability of the oc- 
curence of a "centauro", P(n = 0) = $(0), given by (|49|), is plotted versus 
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the average multiplicity of the system considered for different values of the 
parameter RA, which determines the maximal eigenvalue Ao through the 
relation (fE^). One sees that for RA < 1 this probability remains substan- 
tial even for rather large values of the total multiplicity. With increasing 
-RA, however, it drops pretty fast even at moderate values of the average 
multiplicity. 

6 Discussion and outlook 

One sees from the results of the previous sections that the effects of quantum 
interference can modify substantially the multiplicity distributions expected 
naively from an uncorrelated emission of "distinguishable" particles. The 
modifications become spectacular when the system of particles approaches 
the critical point: the multiplicity distribution becomes very broad and does 
not resemble in any way the original Poisson distribution characteristic of 
uncorrelated emission. This means in particular that one expects relatively 
large probabilities for unusual configurations such as "centauro" or "anticen- 
tauro" events. 

It is thus interesting to discuss the physical conditions for these phenom- 
ena to occur. Considering the example of the Gaussian density matrix, it is 
seen from Fig. 1 that the behaviour of the system is mainly determined by the 
parameter RA and that spectacular effects in charge ratios occur when RA 
is of order 1 or smaller (as already noted, c.f. (p9|), RA > |). The condition 

RA ~ 1 (54) 

implies, for large multiplicity, either a very large particle and energy den- 
sity (when R is small and A takes its "canonical" value of ~ lfm^ 1 ), or 
a "canonical" energy density of about 1 Gev/fm 3 and a very small average 
momentum of the particles in the cm. of the system. To be more specific, a 
system of 100 pions satisfying (|>|) at A = 200 MeV would correspond to the 
energy density of about 35GeV/4fm 3 m 10GeV/fm 3 . On the other hand, for 
an energy density of about 1 GeV/fm 3 (and thus R correspondingly larger), 
A should not exceed 100 MeV. 

Clearly, the probability of creating a " centauro" is enhanced if both effects 
cooperate. We thus conclude that (a) the probabilty of creating a "centauro" 
is enhanced in an environment of high energy density and that the pions 
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emerging from the "centauros" are likely to exibit abnormally small relative 
momenta. 

Finally, let us stress the crucial role played in this analysis by the "orig- 
inal" Poisson multiplicity distribution Po(n) and by the reference density 
matrix po(q,q') for the distinguishable particles. This distribution is easily 
identified when working with Feynman diagrams^, but it is undefined experi- 
mentally. In order to extract the HBT effect from experiment, one would like 
to compare the data with a reference distribution where the HBT effect has 
been switched off. The problem, how to define operationally this reference 
distribution has been discussed for many years without a generally accepted 
conclusion [|l(| |TTJ, for a review c.f. ||. Consequently, it is not possible 
to check separately the assumptions about po(q,q'), Po(n) and the analysis 
of the HBT effect. In the present paper we have chosen simple, but rather 
general, assumptions about po, Pq and concentrated on the HBT effect. No 
doubt, however, more work remains to be done in order to find a realistic 
"distribution of distinguishable particles" po(q, q'), Po{n). 
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Figure caption 

Fig.l. Frequency of "centauro" events as function of the average total 
multiplicity for various values of the parameter RA. 
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